Divergence of the Chaotic Layer Width and Strong Acceleration of the Spatial 
Chaotic Transport in Periodic Systems Driven by an Adiabatic ac Force 
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We show for the first time that a weak perturbation in a Hamiltonian system may lead to an 
arbitrarily wide chaotic layer and fast chaotic transport. This generic effect occurs in any spatially 
periodic Hamiltonian system subject to a sufficiently slow ac force. We explain it and develop an 
explicit theory for the layer width, verified in simulations. Chaotic spatial transport as well as 
applications to the diffusion of particles on surfaces, threshold devices and others are discussed. 
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The basic chaotic formation in perturbed Hamiltonian 
systems is 0, 0, a chaotic layer associated with a 
separatrix of the unperturbed (integrable) Hamiltonian 
system. Even in the simplest case, when the unperturbed 
Hamiltonian system is one-dimensional while the pertur- 
bation is time-periodic, both the transport within the 
layer JH, 0, |(| and its structure on the Poincare section 
l&UiLii relating to the homoclinic tangle, are very com- 
plicated. At the same time, the boundaries of the layer 
are well defined Q as the last invariant curves which 
limit the layer from above and below in the energy scale, 
and may be accurately found in numerical simulations 
For both theory and applications, one of the most 
important characteristics of the layer is its width in en- 
ergy Ul □ III la 111 or in related quantities It 
might be assumed that, if the perturbation is very weak, 
then the layer should necessarily be narrow. This natural 
assumption seems to be su ppo rted by numerous examples 
(e.g. [1SSS0 HUES Ull) However, we show in the 
present Letter that the situation may drastically differ 
(Fig. 1) in a rather general case, namely for any spatially 
periodic system |l^ driven by a slow ac force: the upper 
energy boundary of the layer diverges as the frequency 
of the force goes to zero. We explain this, develop a the- 
ory and verify it by simulations. We also demonstrate 
that the chaotic transport in space may be very fast in 
the adiabatic case, on sufficiently long time-scales, and 
discuss some applications. 

Before moving on to the detailed consideration, we 
comment on the relation to [lo|. It is shown in ^(j that, 
for a system with adiabatically slowly pulsating parame- 
ters, the homoclinic tangle covers most of the range swept 
by the instantaneous separatrix. If the pulsation of pa- 
rameters in [inj | were weak, the range swept would be 
narrow. Our system essentially differs from that consid- 
ered in llOl f our perturbation is not parametric) so the 
result of (TfJ cannot be directly applied to it. If nev- 
ertheless the result of 0] were formally generalized to 
our system, it would give that the homoclinic tangle in 
the adiabatic limit covered the whole phase space, thus 
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FIG. 1: (a). The trajectory in the stroboscopic (for t = 
n2n/ujf with n = 0, 1,2, ...) Poincare section for the system 
(1) with u>o = 1, h = 0.01, ujf = 0.01 and initial conditions 
(p(0) = 0, g(0) = tt). Number of points is 20000. For the sake 
of compactness, we map all points onto the interval [— 7r,7r[: 
q — > q — 2n[(q + n)/(2n)]. This mapping does not affect the 
energy E = p 2 /2 + Uo(q) and, hence, neither the chaotic layer 
width, (b). Spectral dependence of the maximum E max of 
the energy E in the chaotic layer (u>o = 1, h = 0.01). 



hinting at a divergence of the chaotic layer width. 

Periodic systems and ac forces are widespread in na- 
ture. An archetypal example is a pendulum driven by 
a weak single-harmonic ac force: 



Q = P, P 
U = U (q) = 
U = U(q,t)-- 



- —dllo/dq — Iiujq sm(ujft) 
-Wo cos(q), 

U (q) + qhwl sin(cj/i), 



-dU/dq, 



(1) 



/i < 1. 



To get some insight into the physical origin of the phe- 
nomenon, and to obtain a qualitative explanation of Fig. 
1, we first consider the strong adiabatic limit 



Ljf/(hu Q ) < 1, 



(2) 



a much wider range of Wf will be considered afterwards. 

Given that uif is small, the system may be considered 
as moving in the quasi-stationary "wash-board" poten- 
tial U with slope Iiluq sm(uj ft) . The absolute value of 
the slope is ~ Iiluq for most of the time while its sign 
is positive during odd half-periods of the perturbation, 
changing to the opposite one in the even half-periods. If 
the system is initially at the top of the potential bar- 
rier, then, owing to the condition (2), even a small slope 
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~ Jiu)q is sufficient to accelerate it during the first half- 
period up to a negative velocity of large absolute value. 
In the second half-period, the slope changes its sign, re- 
sulting in a braking effect, so that the velocity returns 
close to zero at the end of the period. Assuming that the 
maximum of the kinetic energy K, = p 2 /2 greatly exceeds 
the spatial modulation of U(q,t), 

K-max 3> w , (3) 

one may neglect the term Uo(q) in U(q, t) while describing 
the major part of the trajectory (namely, when IC S> u>q) 
during the first period of perturbation: the equations of 
motion reduce to those of a free particle driven by the 
time-periodic force, which are solved exactly: 

q(t) = A+(b- ^] t + ^ sinfw/t), 
V Uf J ujj 

hi ? 

p(f) = g = fl-^(l-cos( W/ t)), (4) 

where A = q(Q) and B = p(0). 

For the given initial state (p(0) = 0, q(0) = tt), B = 
so that the velocity (4) oscillates from to —2hu>Q/u>f 
while the kinetic energy JC = p 2 j2 oscillates from to 

K. m ax = 2h 2 uj£/uj 2 f = 2uj 2 /e 2 , e < 1. (5) 

Eq. (5) confirms the validity of the assumption (3), and 
the maximum of the energy does diverge as Uf — > 0. 

But Eqs. (4), (5) do not explain chaos in the system. 
Like in other cases 0, 0, S 0], the chaos onset in our 
system is related to the motion near the unperturbed 
separatrix. Its rigorous treatment is complicated (cf. the 
conventional adiabatic case |lo||l but is not essential for 
the quantity of main interest in our paper, i.e. for the 
chaotic layer width: the width is determined mainly by 
the "regular" parts of the chaotic trajectory (described 
by (4)). So, we describe the chaos onset just qualitatively. 

Define h n = h n (u>f) as the value of h for which the sys- 
tem starting from the top of the barrier (p(0) = 0, q(0) — 
tt) arrives in the end of the first period of the perturbation 
t^ 1 ' = 2tt /u)f at the top of the nth barrier i.e. q(t^) = 0, 
q(t^) = tt — 2im where n is a large positive integer. If 
h = h n , the velocity of the system at the instant when 
it passed the previous (i.e. (n — l)st) barrier top can be 
shown to be — q c with q c ~ ^hujfUjQ ln(wo / '(hu> /)) <C Wo- 

Consider h slightly smaller than h n . For any t, —q(t) 
on the regular part of the trajectory is slightly smaller 
than that for h = h n . Hence, the velocity on the true 
trajectory becomes zero slightly to the right of the nth 
barrier top slightly before i.e. the trajectory reflects 
from the right slope of the nth barrier and qft 1 ') > 0. 
The new round of acceleration to the left starts only after 
the next reflection (from the left slope of the (n — l)st 
barrier) i.e. with a small delay At\ with respect to fW; 



Ah ~ uo Q 1 \n{ul/ql) ~ uo Q 1 In(w /(ftw/)) < u f 1 . (6) 

This new round of acceleration-braking is described by 
Eq. (4) with B determined by the condition p(2ir /ojf + 
Aii) = 0, i.e. B = (w /e)(l - cos(w/Aii)) > 0. 

Consider h slightly larger than h n . For any t, —q(t) is 
slightly larger than that for h = h n on the regular part of 
the trajectory. Hence, at t^\ the system arrives slightly 
to the left of the top of the nth barrier while moving to 
the left and with an energy close to the top of the barrier 
level. When the coordinate of the system approaches the 
vicinity of the top of the next (i.e. (n + l)st) barrier of 
U(q, t), the latter becomes sufficiently lower than the top 
of the nth barrier and the system passes over it rather 
than reflects. Thus, the second round of acceleration- 
braking is described by Eq. (4) with B = q(t^') being 
negative with a small absolute value. 

At the end of the second acceleration-braking round, 
the system may again either pass over a barrier or re- 
flect from it and even pass in the backward direction 
for a few periods of the potential. In the latter case, 
the new acceleration-braking round is additionally de- 
layed with respect to the perturbation by some time Ai2 : 
this round is described by (4) with yet a larger posi- 
tive B = (u> /e)(l - cos(w/[A*i + Ai 2 ])) > 0. As time 
goes on, the process of reflection develops and, on long 
time-scales, the delay due to the reflections is accumu- 
lated. So, B gradually growes in average until it gets 
close to B max = 2wo/e, then gradually decreases in av- 
erage until it gets close to 0, etc. The sequence of reflec- 
tions/passings is random unless h = h n , but even in this 
case the sequence is random if the initial state is shifted 
from the top of the barrier. As reflections/passings and 
intervals Aij are random, variations of B are random too: 
so, the trajectory is chaotic on large time-scales. 

Note also that the action I = (2n)^ 1 § p dq 16], con- 
ventionally 0, 0, Q] chosen as the lowest-order adia- 
batic invariant, is not conserved for motion above the 
barrier in our system: on the major part of a trajectory, 
K 3> uj 2 and hence / « p while \p\ varies in a wide range 
(from to ~ wo/e, for the chaotic trajectory). The cor- 
rect lowest-order adiabatic invariant for our system is 

hi > 2 

/ = /+-=^(l-C08(«/t)). (7) 

w / 

To the lowest order in e, I coincides with the integration 
constant B of (4) on the major part of the trajectory [l7| . 
So, the chaotic layer width in / is the same as in B, i.e. 
equal to B max = 2htjQ/u)f, diverging as ujf — > 0. 

Let us consider the problem of the layer width in a 
more general case, when Eq. (2) may not be satisfied. As 
explained above, the trajectory that can both pass over 
a barrier and reflect from it is chaotic. But the strong 
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inequality (3) for IC max may not hold now (cf. Eq. (5)) 
so the term Uo in U cannot be neglected even on regular 
parts of the chaotic trajectory. Let the adiabatic approx- 
imation be still valid (the explicit condition is derived in 
Eq. (18)). As h is small, the turning point of the chaotic 
trajectory is situated near the top of a potential barrier: 



p(t s ) = 0, q(t s ) w 7r - 27rm, 



(8) 



where m is integer and t s is one of instants when p = 0. 

Let IC n be the kinetic energy K. at the instant t n when 
the system crosses the coordinate of the top of the nth 
barrier of the auxiliary potential U, i.e. 



7r — 27m, 



ICr. 



P 2 (t n ) 



(9) 



The coordinate of the next barrier top crossed by the 
system is 



Qn+5„ 



2nSr, 



- S iga(p(t n )). (10) 



In the adiabatic approximation, the change of K. during 
the time t n+ s n — t n is the following: 



K 



n+8 n 



fCr, 



-{o.n+S„ ~ q n )hu>Q sm(ujft n ). (11) 



The interval of time t n+ g n — t n can be evaluated in the 
adiabatic approximation as follows 



t n +s n - 1„ — — (12) 

f q " dq 2ViK(x) 

A„ +i „ 6 n y/2(JC n + + cos(q))) 

K(at) = 



dip 



1 



1 — x sin (ip) 



l + fc»/(2w§)' 



K(a;) is the full elliptic integral of the 1st order 

Significant changes of /C n take place on the time scale 
luJ 1 while the adiabatic condition means that 

tn+5 n — t n <C ^J 1 - (13) 

So, discrete changes of K, n = K,(t n ) may be replaced by 
the continuous change of a function K,{t n ) with derivative 

d/C ^ fcn+8 n — K-n 
dt n 



■KhujlJl + fc/{2ojl) 
w -sign(p(t n )) - sm{u f t n ). (14) 

K(l/(l + /C/(2^))) 

Separating variables, evaluating the resulting integrals 
and taking into account that, within the present simpli- 
fied description of motion, sign(p) changes when fC = 0, 
we obtain the transcendental equation for K,: 



E(x) nh loq 

—f=- = 1 + —. cos(w/t s ) - cos(w/t n ) , (15) 



E{x) = J dip ^/l-zsin (ip), 2r=ll + ^-|J . 

Here, E(x) is a full elliptic integral of the 2nd order (lif . 

Due to the random-like changes of t s , cos(uift s ) on long 
time-scales is densely distributed over the range [—1,1]. 
Maximizing K, with respect to t s and t n and taking into 
account that E(t n ) = JC(t n ) + Uo(q(t n )) w K.(t n )+ojQ, we 
finally obtain the transcendental equation for the upper 
energy boundary of the chaotic layer E max : 



E(x) _ 7rfe loq 

y/x 2 ' 



2wg 



(16) 



E m ax monotonously decreases from oo to Wq as u)f/hu>o = 
e increases from to oo. Fig. 2 shows that Eq. (16) nicely 
describes the simulations in a wide range of ujf. 




FIG. 2: Spectral dependence of the maximum energy in the 
chaotic layer (uiq = 1, h = 0.01): circles, solid, dashed and 
dotted lines correspond respectively to simulations, the gen- 
eral formula (16), the asymptotes (5) and (17). 

For e <C 1, the root (E max — lUq)/lUq of Eq. (16) is 
large. So, x w 2u)%/E max — > while E(x — > 0) — > tt/2, 
and the solution of (16) reduces to the asymptote (5). 



For e ^> 1, the root (E r , 



loD/loI of Eq. (16) is 



small. Using the asymptote 18] for E(x — > 1), Eq. (16) 
can be reduced to the following asymptote for E max : 



E n 



32 



aln(et) 



— , e>l. 

7T 



(17) 



Using the estimate (12) and the asymptote K(x — > 
1) = 0.51n(16/(l-cc)) [Til and allowing for Eq. (17), we 
may express the adiabatic condition (13) explicitly: 



uj f < u /hi(l/h). 



(18) 



Let us discuss chaotic transport. For most physical 
applications, transport in coordinate is relevant. For 
sufficiently long time-scales, the chaotic trajectory pro- 
vides large-scale displacements in both directions, un- 
like regular trajectories. Generally, the mean-square dis- 
placement for chaotic transport depends on time as 
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((?(*) ~ <?(0)) 2 ) = Dt b with < b < 2 ((...) means av- 
eraging of the initial conditions over the chaotic layer). 
The larger D and b are, the faster the transport is. It is 
suggested in 0], basing on numerical results, that b — > 1 
as LOj — > 0. As for D, we suggest that it strongly diverges. 
Indeed, chaotic trajectories generally spend most of the 
time close to the boundaries of regions of regular motion 
6J. In our case, it means that, for most of the time, the 
trajectory moves close to either the upper border of the 
chaotic layer, with the average velocity v w hur^juif, or 
the lower border, with the average velocity —v. We call 
such regimes acceleration-braking flights (cf. Levy flights 
0,H3)i distinguishing them from the regime of diffusion 
across the layer. The duration tf of the flight may be es- 
timated from the analysis of the diffusion of B near the 
boundary of the layer. The diffusion constant for B may 
be roughly estimated as Db ~ ((A£?) 2 )/(27r/w/) where 
((A£>) 2 ) is the average squared change of B at the end 
of a driving period: ((AS) 2 ) ~ B 2 max (uj /Ati) 4 . Then, 



Dj 



In (u) /(hu)f)) 



(19) 



Finally, D may be estimated as the ratio between the 
squared length of the flight / 2 and its duration tf. 

(uJo/ujf) 7 



D ~ -i- 



tf 



v 2 t f 



Ulgh 



In (u> /(hu)f)) 



(20) 



The above analysis provides intuitive arguments in favor 
of a strong acceleration of the chaotic transport in space 
as ujf — > 0, and simulations support this. Still, a thor- 
ough numerical study and a rigorous evaluation of D(ujf) 
as well as a proof that b(uif — > 0) = 1 are necessary. 

Apart from purely dynamic phenomena, our work may 
have a strong impact on noise-induced phenomena: e.g. 
diffusion of particles on surfaces (see [ijj and references 
therein), that plays an important role in many mod- 
ern technologies involving self-assembled molecular-film 
growth, catalysis, and surface-bound nanostructures |2lj . 
Numerous studies of atoms (e.g. |22j), or gan ic molecules 
(e.g. |23(), and even metal clusters (e.g. 24]) show that 
the long jumps may play a dominant role in the diffu- 
sion: this means that the damping is small. If we add a 
weak ac drive, then transient chaos Q arises in a region 
of phase space which approximately coincides with the 
chaotic layer developed for zero damping (cf. Q). In 
the light of the present results, the latter means: if the 
driving is adiabatic, then, as soon as the noise activates 
the particle to the lower boundary of the layer (i.e. ap- 
proximately to the barrier level), the further transport 
is provided by fast (compared to slow noise-induced dif- 
fusion) chaotic transport. Thus, even a weak drive may 
drastically increase the speed of the diffusion. 

One more possible application concerns a threshold de- 
vice (cf. [13 )• If the device was based on a noise-driven 



spatially periodic system and switched at the energy level 
E t h significantly exceeding the barrier level Ub, then the 
addition of even a weak but sufficiently slow ac drive 
would give rise to a marked decrease of the activation 
energy - from E t h to Ub, thus leading to a drastic in- 
crease of the flux. This may also be used as a sensitive 
method to measure the amplitude/frequency of the drive. 

In conclusion, we have discovered the divergence of the 
chaotic layer width and the related strong acceleration 
of the chaotic spatial transport in ac driven periodic sys- 
tems in the adiabatic limit. The mechanism is a combina- 
tion of acceleration-braking flights of high average kinetic 
energy and of small random changes of the adiabatic in- 
variant near the separatrix. Applications to diffusion on 
surface and to threshold devices have been suggested. 
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